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Lightning flashes produce nitric oxide (NO) which is an important precursor for ozone in 
the lower atmosphere. However, the amount of production of NO by lightning globally is very 
uncertain, and in fact has the largest uncertainty of all of the sources of NO both man-made and 
natural. Better estimates of the lightning NO production are needed such that improved 
estimates of the role of natural and man-made sources of NO in the chemical production of 
ozone in the atmosphere can be made. Ozone, at the altitudes where lightning NO is an 
important contributor, is an effective absorber of infrared radiation, and therefore climatically 
important. Since satellite measurements now provide us with good estimates of the number of 
lightning flashes occurring globally, most of this uncertainty stems from the large range in 
estimates of the NO production per flash. In this paper we use a combination of models and 
observations to deduce the average NO production per cloud-to-ground (CG) flash and per 
intracloud (IC) flash during a thunderstorm that was studied by aircraft, radar, and lightning 
sensors over southern Germany in 1998. The storm was simulated with a cloud model, and the 
cloud model data were used to run a cloud-chemistry model that includes a lightning NO 
generation scheme. Observed lightning flash rates are input to the model and the model places 
each flash into appropriate regions of the storm. The winds in the model carry the NO and other 
trace gases through the cloud, while chemical reactions involving NO, ozone, and other species 
continuously occur. Various combinations of production of NO by CG and IC flashes are tested 
in the model and the results of each combination are compared with observations of NO taken by 
aircraft in the anvil of the thunderstorm. The NO production scenario with the best match to the 
observed NO concentrations is considered to be the most likely. For this storm a scenario of NO 
production with the average production per CG flash equal to that of an IC flash yielded the best 
results. This outcome along with similar results for a US Great Plains storm are in contrast to 
older literature which proposed that an IC flash makes only 10% as much NO as a CG flash. 
However, the IC and CG NO production in the analyzed storm were both less than the amount 
assumed for a CG flash based on the literature. These results have implications for the vertical 
distribution of the lightning NO production since much of the IC lightning activity occurs in the 
middle and upper portions of the cloud. Ozone within the storm cloud decreased during the 
simulation due to reactions with the very strong NO plume generated by lightning. However, 
downwind of the storm, ozone production occurs as the NO plume dilutes. 
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ABSTRACT 


1 
2 

3 The July 21, 1998 thunderstorm observed during the European Lightning 

4 Nitrogen Oxides Project (EULINOX) project was simulated using the three-dimensional 

5 Goddard Cumulus Ensemble (GCE) model. The simulation successfully reproduced a 

6 number of observed storm features including the splitting of the original cell into a 

7 southern cell which developed supercell characteristics, and a northern cell which became 

8 multicellular. Output from the GCE simulation was used to drive an offline cloud-scale 

9 chemical transport model which calculates tracer transport and includes a 

10 parameterization of lightning NO x production which uses observed flash rates as input. 

1 1 Estimates of lightning NO x production were deduced by assuming various values of 

12 production per intracloud and production per cloud-to-ground flash and comparing the 

13 results with in-cloud aircraft observations. The assumption that both types of flashes 

14 produce 360 moles of NO per flash on average compared most favorably with column 

15 mass and probability distribution functions calculated from observations. This assumed 

16 production per flash corresponds to a global annual lightning NO x source of 7 Tg N yr' 1 . 

17 Chemical reactions were included in the model to evaluate the impact of lightning NO x 

18 on ozone. During the storm, the inclusion of lightning NO x in the model results in a 

19 small loss of ozone (on average less than 4 ppbv) at all model levels. Simulations of the 

20 chemical environment in the 24 hours following the storm show on average a small 

21 increase in the net production of ozone at most levels resulting from lightning NO x , 

22 maximizing at approximately 5 ppbv day' 1 at 5.5 km. Between 8 and 10.5 km, lightning 

23 NO x causes decreased net ozone production. 
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1 1. INTRODUCTION 

2 The nitrogen oxides NO and NO 2 (NO x ) are important trace gases in the 

3 troposphere, the presence of which in sufficient quantities may result in the 

4 photochemical production of ozone. Fossil fuel combustion, biomass burning, microbial 

5 activity in soils, and lightning are considered the four major sources of tropospheric NO x 

6 [Bradshaw et al., 2000]. Though lightning produces much less NO x than the 

7 anthropogenic sources, it is a particularly significant source because lightning produces 

8 NO x mainly in the middle and upper troposphere where it is longer lived, typically at 

9 more dilute concentrations, and, consequently, more efficient at producing ozone than in 

10 the boundary layer where the majority of NO x is emitted. 

1 1 The strength of the global NO x source from lightning remains uncertain, with 

12 values in the literature ranging over an order of magnitude (3.2-26 Tg N yr' 1 ) [Bradshaw 

13 et al., 2000], although most recent estimates are confined to the lower half of this range 

14 as summarized in Zhang et al. [2003]. Contributing to the widely varying global 

15 estimates are uncertainties in both the total number of flashes and amount of NO x 

16 produced per flash or per meter of flash length. There is also debate regarding the 

17 relative NO production by intracloud (IC) and cloud-to-ground (CG) flashes (e.g. Ridley 

18 et al., 2005). On the basis of previous studies which suggested that IC flashes were less 

19 energetic than CG flashes (e.g. Holmes et al., 1971), many studies of lightning NO x 

20 production have assumed that production by an IC flash (Pic) is less than production by a 

21 CG flash (Pcg)- Price et al. [1997] assumed that production by an IC flash (Pic) was one 

22 tenth of production by a CG flash (Pcg) in calculating global lightning NO x production. 

23 Gallardo and Cooray [1996] suggested that IC flashes may dissipate nearly as much 

24 energy as CG flashes and therefore Pic may be on the order of Pcg- Supporting the 
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1 Gallardo and Cooray [1996] hypothesis, a 2-D cloud-scale modeling study by DeCaria et 

2 al. [2000] suggested that the Pic/Pcg ratio is likely between 0.5 and 1.0, and a 3-D 

3 simulation of the same storm narrowed this range to between 0.75 and 1.0 [DeCaria et al., 

4 2005]. 

5 Estimates of NO x production by lightning using laboratory experiments, 

6 theoretical assumptions regarding the physics of lightning flashes, data from field 

7 experiments, and global models were summarized in Zhang et al. [2003] and range from 

8 8 to 5000 moles of NO per flash. Below we present estimates of NO x production per 

9 flash and per meter flash length obtained by various investigators using data from several 

10 recent field projects including STERAO (Stratosphere Troposphere Experiment: 

1 1 Radiation, Aerosol, Ozone) over northeastern Colorado in 1996, LINOX (Lightning- 

12 produced NO x ), in southern Germany in 1996, and EULINOX (The European Lightning 

13 Nitrogen Oxides Project) over central Europe in 1998. Estimates of NO x per lightning 

14 flash are summarized in Table 1 while estimates of NO x production per meter flash length 

15 are presented in Table 2. For the July 12 STERAO storm over Colorado, DeCaria et al. 

16 [2000] used a 2-D cloud-resolving model, and DeCaria et al. [2005] used a 3-D cloud- 

17 resolving model, anvil NO observations, and observed lightning flash rates calculated 

18 from interferometer and NLDN (National Lightning Detection Network) observations to 

19 estimate Pcg and the likely Pic/Pcg ratio. Based on analysis of lightning interferometer 

20 and aircraft data, Skamarock et al. [2003] estimated NO x production per interferometer 

21 flash and per meter flash length produced by lightning in the July 10 STERAO storm. 

22 The interferometer flash data included many short duration flashes that likely would not 

23 be detected by other lightning detection systems. Whether or not these short duration 

24 flashes are productive of NO is an open question. If not, then the NO production per 
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1 meter in the July 10 STERAO storm would be larger than computed by Skamarock et al. 

2 [2003]. Stith et al. [1999] estimated NO production per meter flash length based on an 

3 analysis of NO spikes measured by aircraft during the same storm, as well as other 

4 STERAO storms. Using a similar approach, Huntrieser et al. [2002] estimated NO 

5 production per meter flash length for the July 21 EULINOX storm over Germany. For 

6 the same storm, an average production of NO per flash was estimated by Thery et al. 

7 [2000] based on an average flash length of approximately 30 km calculated from 

8 interferometer data (mostly IC flashes). Through analysis of CG lightning peak current 

9 data for the July 21 EULINOX storm as recorded by BLIDS (Blitz Informationsdienst 

10 von Siemens), Fehr et al. [2004] estimated Pcg and a Pic/Pcg ratio of 1.4 which was 

1 1 confirmed by a cloud-resolving model simulation that included a lightning 

12 parameterization, followed by comparison with anvil NO x observations. Based on 

13 aircraft chemical measurements from LINOX, Huntrieser et al. [1998] estimated NO 

14 production per flash while Holler et al. [1999] estimated the amount of NO produced per 

15 meter flash length for the same project using lightning observations from LPATS 

16 (Lightning Position And Tracking System). All of the above estimates of NO production 

17 are based on aircraft NO x measurements in storm anvils. In contrast, Langford et al. 

18 [2004] used remote sensing techniques to measure NO 2 column abundances in the core 

19 region of a thunderstorm over Boulder, Colorado, and estimated a larger value of 963 

20 moles of NO produced per CG flash. The wide range of values found using data from 

21 various field projects suggests that lightning NO x production may vary greatly between 

22 storms and individual flashes, though all estimates of lightning NO x production contain 

23 significant uncertainty because of the number of assumptions required in translating 

24 aircraft observations into NO production values. Further investigation is required to 
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1 constrain the magnitude of NO production per flash in order to reduce uncertainty in the 

2 global source estimate. 

3 In this study, the three-dimensional (3-D) Goddard Cumulus Ensemble (GCE) 

4 model is used to simulate the dynamic evolution of the EULINOX thunderstorm 

5 observed on 21 July 1998. This output is used to drive an offline cloud-scale chemical 

6 transport model (CSCTM) which includes a parameterized lightning NO x source. The 

7 CSCTM results are first used to estimate NO x production per CG flash and per IC flash. 

8 The lightning parameterization used in this model represents an advancement over the 

9 bulk approach that was developed by DeCaria et al. [2005], in that the NO production of 

10 individual flashes is simulated in the new scheme. In an alternative parameterization, 

1 1 interferometer estimates of average hourly flash length during the storm are employed in 

12 the model to yield an estimate of NO production per meter flash length. The ability to 

13 estimate NO production per unit flash channel length also represents an advancement in 

14 the model over the version used in earlier work (e.g. DeCaria et al., 2005). Chemical 

15 fields from the CSCTM are used in a chemistry-only version of the model to estimate 

16 downstream ozone production in the 24 hours following the storm. 

17 Fehr et al. [2004] also simulated the July 21 EULINOX storm using a cloud 

18 resolving model with a parameterized lightning NO x source to investigate lightning NO x 

19 production and transport, but did not examine the transport of other species, or the 

20 chemical impact of lightning NO x during or after the storm. Lightning NO x emissions 

21 were represented by a Lagrangian particle model with particles distributed within a 

22 vertical lightning channel which does not account for the tortuosity of an actual lightning 

23 flash. The Fehr et al. [2004] study also did not attempt to quantify production per meter 

24 of flash channel length. In contrast, the research reported here considers the transport of 
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1 other species, the production of NO x per meter of flash channel length, and the resulting 

2 photochemical ozone production/loss. 

3 Section 2 of this paper provides background information on the EULINOX 

4 project and the investigated thunderstorm, as well as the models used. Section 3 

5 discusses the results of the various simulations, and section 4 contains a discussion of 

6 conclusions drawn from this research. 

7 2. OBSERVATIONS AND NUMERICAL MODELS 

8 2.1. Overview of the EULINOX project and the July 21, 1998 storm 

9 The European Lightning Nitrogen Oxides Project (EULINOX; Holler and 

10 Schumann, 2000; Huntrieser et al., 2002) was conducted in central Europe during June 

1 1 and July 1998 with the goal of better understanding lightning NO x production. During 

12 the project, airborne measurements were collected by the Deutsches Zentrum fur Luft- 

13 und Raumfahrt (DLR) Falcon and Do228 research aircraft. Both measured NO x , O3, CO 

14 and CO2, as well as meteorological parameters in and around thunderstorms in the region 

15 of the experiment. The Do228 flew primarily in the boundary layer and lower 

16 troposphere below 4 km, while the Falcon investigated the upper troposphere and 

17 performed a number of anvil penetrations through monitored thunderstorms. CG 

18 lightning occurrences were recorded by an LPATS system known as BLIDS. In addition, 

19 total 3-D lightning activity (IC+CG) within the 100 km area surrounding the EULINOX 

20 operation center was mapped by a VHF (very high frequency) interferometer from the 

21 French Office Nationale d’Etudes et de Recherches Aerospatiales (ONER A). Both radar 

22 and satellite observations were used to monitor the development of thunderstorms in the 

23 region. 
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1 On the evening of 21 July 1998 the evolution of a severe thunderstorm west of 

2 Munich, Germany was documented. Figure 1 [Holler et al., 2000] shows the evolution of 

3 the thunderstorm as observed by the Polarisation Diversity Radar (POLDIRAD) radar at 

4 1° elevation. After an initial period of intensification, the storm split into two distinct 

5 cells evident on the radar image observed at 1852 LT (1652 UTC). The northernmost 

6 cell became multicellular in structure and was observed to decay soon after the cell- 

7 splitting event, while the southern cell strengthened and developed supercell 

8 characteristics including a distinct hook echo (see radar image observed at 1933 LT; 

9 Holler etal., 2000). 

10 From 1400 to 2200 UTC total lightning activity within the entire storm was 

1 1 monitored by the ONERA interferometer, and cloud-to-ground lightning activity was 

12 recorded by BLIDS. Thery [2000] compared flashes recorded by BLEDS and the 

13 ONERA interferometer and found that many low amplitude positive flashes recorded by 

14 BLIDS were IC rather than CG. This was also true for flashes with weak negative peak 

15 current, but to a lesser degree. Flashes with peak current between -5 and +15 kA are 

16 removed [Holler et al., 2000] and the remaining BLIDS flashes counted to obtain CG 

17 flash rates, which were subtracted from counts of interferometer flashes to yield IC flash 

18 rates. By comparing plots of the locations of interferometer and BLIDS flashes with 

19 plots of radar reflectivity, flashes belonging to the northern and southern cells were 

20 identified. Figure 2 shows the estimated flash rates for the northern and southern cells 

21 following the cell splitting. The southern cell contained the majority of lightning activity 

22 with IC flashes dominating the total lightning activity after 1710 UTC. In the period 

23 shown, 360 CG flashes and 2565 IC flashes were recorded in the southern cell (mean 
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1 IC/CG ratio = 7.1) while in the northern cell there were 289 CG flashes and 815 IC 

2 flashes (mean IC/CG ratio = 2.8). 

3 The storm was penetrated seven times during the period from 1735 to 1842 UTC 

4 by the Falcon while flying between 6.3 and 9.2 km AGL near the active convective cells 

5 and in the lower anvil. Figure 3 [Huntrieser et al., 2002] shows the Falcon flight track 

6 superimposed on a map of VHF signals recorded by the interferometer from 1740 to 1810 

7 UTC. During this period, nearly 400 IC and CG flashes were recorded. The Falcon flew 

8 between 8 and 9 km AGL in the vicinity of the maximum lightning activity where the 

9 majority of freshly produced NO would be observed. NO mixing ratios up to 25 ppbv 

10 were observed during this period [Huntrieser et al., 2002]. 

1 1 Fehr et al. [2004] simulated the EULINOX storm of July 21, 1998 using a 

12 modified, cloud-resolving version of the PennState/NCAR Mesoscale Model 5 (MM5) 

13 which included a parameterized lightning NO x source with emissions represented by a 

14 Lagrangian particle model. Lightning NO x particles were distributed within a vertical 

15 flash channel with the vertical distance between particles depending on atmospheric 

16 pressure. Flash rates were parameterized using the methods of Price and Rind [1992, 

17 1993] and Pickering et al. [1998] and compared with observed flash rates. The 

18 parameterization overestimated total lightning activity considerably, so observed flash 

19 rates were used to estimate lightning NO x production by comparing model results with 

20 NO measurements taken during three of the seven Falcon anvil penetrations. A 

21 production scenario in which a CG flash produces approximately 330 moles of NO per 

22 flash and an IC flash is 1.4 times more productive of NO than a CG flash compared 

23 favorably with the observations used. An estimated 50-80% of the lightning produced 
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1 N0 X was transported to the anvil region, with 97% of anvil NO x resulting from IC 

2 flashes. 

3 2.2. Model Description 

4 In order to simulate the July 21 EULINOX storm, the 3-D Goddard Cumulus 

5 Ensemble Model (GCE), a non-hydrostatic, cloud resolving model, was employed with a 

6 horizontal resolution of 2 km and vertical resolution of 0.5 km. The GCE uses a Kessler- 

7 type scheme for cloud water and rain [Kessler, 1969; Houze, 1993], and a three-category 

8 scheme from Lin et al. [1983] for cloud ice, snow, and hail to parameterize cloud 

9 microphysics. The open boundary conditions of Klemp and Wilhelmson [1978] are 

10 employed at the lateral boundaries. A complete description of the GCE model is given in 

1 1 Tao and Simpson [1993], with updates described by Tao et al. [2001]. 

12 An associated offline 3-D Cloud-Scale Chemical Transport Model (CSCTM) was 

13 developed at the University of Maryland [DeCaria et al., 2000; DeCaria et al., 2005], 

14 Temperature, density, wind, hydrometeor (rain, snow, graupel, cloud water, and cloud 

15 ice), and diffusion coefficient fields from the 3-D GCE are read into the model every ten 

16 minutes in the simulation, and these fields are then interpolated to the model time step of 

17 15 seconds. The transport of chemical tracers is calculated using the van Leer advection 

18 scheme. The CSCTM employs parameterizations of lightning NO x production which use 

19 observed flash rates, avoiding the difficulty of explicit modeling of thunderstorm 

20 electrification. A passive version of the CSCTM includes only the transport of tracer 

21 species and production of lightning NO x without any chemical reactions. In the passive 

22 version, the IC- and CG- produced NO x can be isolated from pre-existing NO x which 

23 allows a number of different production scenarios to be evaluated quickly. To account 

24 for the effects of chemical reactions, the full version of the CSCTM combines tracer 
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1 transport and lightning production with a chemical solver and photochemical mechanism, 

2 better simulating the actual chemical environment within the storm. The same chemistry 

3 reaction scheme as used by DeCaria et al. [2005] is employed here, except that reaction 

4 schemes for isoprene and propene were added. Soluble species are removed from the gas 

5 phase by cloud and rain water. Multiphase reactions are not included, and photolysis 

6 rates are calculated as a function of time. A chemistry-only version of the CSCTM is 

7 used to estimate 24-hour ozone production in the convective outflow that is assumed to 

8 be transported downwind. In this version, chemical reactions and diffusion are included, 

9 while cloud-scale advection is turned off. Clear sky photolysis rates are assumed and the 

10 photolysis rates are calculated as a function of time. 

1 1 Two different parameterizations are used in the model to estimate the production 

12 of lightning NO x . These parameterization schemes differ from the lightning NO x 

13 parameterization used in DeCaria et al. [2005] in which lightning NO x was distributed 

14 bimodally in the vertical and uniformly to all grid cells within the 20 dBZ contour of the 

15 cloud at each level as if the NO x was instantly diffused throughout this region of the 

16 cloud. The DeCaria et al. [2005] model results were compared to the general profile 

17 shape and integrated column mass of observed NOx. This approach performed well for 

18 the July 12 STERAO storm in which the anvil observations were located relatively 

19 distant from the convective cores. In that case, the aircraft measured the integrated 

20 effects of many flashes on NO x mixing ratios. However, in the EULINOX storm the 

21 aircraft flew in a much more electrically active part of the storm, necessitating a different 

22 approach. These new parameterizations attempt to more realistically replicate actual 

23 flashes and the range of NO x mixing ratios observed by putting lightning NO x from 

24 individual flashes into smaller subsets of grid cells within the cloud. 
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1 The first parameterization allows an estimate to be made of NO production per 

2 flash and the second allows an estimate of production per meter flash length. In the first 

3 scheme, observed flash rates for the northern and southern cells are input along with a 

4 scenario of IC and CG production specified in terms of moles of NO produced per flash. 

5 The average horizontal extent of a flash is calculated from interferometer data and is 

6 input for each 3-minute lightning time step, as is the total number of IC and CG flashes in 

7 each cell of the storm. The areas in which lightning occurred in the northern and 

8 southern cells were estimated from plots of observed IC and CG flashes. Areas of 

9 approximately this size were centered 10 km downwind of the maximum updraft of the 

10 northern and southern cells in the model because Holler et al. [2000] noted that in this 

1 1 storm, based on an analysis of interferometer and radar observations, flashes tended to 

12 occur downwind of the updraft. The distance of 10 km was chosen based on visual 

13 inspection of plots of radar reflectivity overlaid with flash locations presented in Holler et 

14 al. [2000]. 

15 The vertical distribution of IC flash channel segments was derived from two 

16 Gaussian distributions, one centered at -30° C and the other at -15° C, which were 

17 summed, while the vertical distribution of CG flash channel segments consisted of a 

18 single Gaussian distribution centered at -15° C [DeCaria et al., 2005]. These 

19 distributions, f(z), determine the number of grid cells in the horizontal to be included in 

20 an IC or CG flash at each model level as shown in Figure 4a. The vertical distribution of 

21 the number of grid cells included in a CG flash results in no grid cells in the lowest two 

22 layers of the model receiving direct placement of lightning NO. To a first approximation, 

23 this configuration is supported by data from 3-D lightning mapping systems which record 

24 a nearly negligible amount of flash channels near the surface compared with that which 
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1 occurs aloft. At the top of the cloud, as determined by the uppermost nonzero value of 

2 f(z), an initiation point is selected at random within the designated area downwind of the 

3 the updraft. After this point is selected, the flash is constrained to an area equaling the 

4 average horizontal flash extent. At each level, the locations of a number of grid cells 

5 given by f(z) are selected at random, such that tortuosity of the flash is simulated. NO 

6 production is distributed to all grid cells along each flash with a dependence on pressure 

7 as described in DeCaria et al. [2000] because of laboratory experiments showing a linear 

8 relationship between pressure and NO production [Wang et al., 1998]. Figure 4b shows a 

9 schematic diagram of lightning NO placement in the southern cell. The dashed line 

10 represents the 20 dBZ contour at 9 km 150 minutes into the simulation. The maximum 

1 1 updraft velocity location is identified by the triangle. The larger box centered 

12 downstream of the updraft velocity maximum designates the area from which an 

13 initiation point for the flash, marked with an open circle, is selected. The smaller box 

14 centered about the initiation point is the area determined by the average horizontal flash 

15 extent. Various NO x production scenarios are simulated to determine which most 

16 closely matches observed NO x mixing ratios in the electrically active region. 

17 The second parameterization scheme is similar to the first, except that the 

18 specified production is per meter of flash channel length and production per flash is 

19 calculated by the model using the average hourly length per flash as given in Thery et al. 

20 [2000]. Lightning flashes are constructed in the same manner as in the first 

21 parameterization, and the most appropriate production per meter flash length is estimated 

22 by comparing results from various production scenarios with aircraft observations. 

23 3. RESULTS 

24 3.1. Cloud simulation 
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1 The GCE model was initialized with a single sounding that included data from a 

2 German Weather Service radiosonde, the ascent of the DLR Falcon aircraft, and a 

3 dropsonde released during the Falcon’s flight, all of which were no more than 90 minutes 

4 ahead of the storm. CAPE for this sounding was 1590 J/kg [Fehr et al., 2004]. The 

5 sounding also shows winds veering at low levels which is conducive to splitting. 

6 Convection was initiated with a warm thermal perturbation. The model domain was 

7 360x328 km, with a horizontal resolution of 2 km. There were 50 vertical levels, with a 

8 resolution of 0.5 km. The model was run for 6 hours. 

9 The GCE simulation successfully reproduced a number of features of the 

10 observed storm. A single cell first appears 20 minutes into the simulation. At 70 

1 1 minutes, the cell splitting begins. Because the early stages of the cell splitting is 

12 observed on radar at ~ 1650 UTC (Figure 1), 70 minutes in the simulation was chosen to 

13 correspond to this time for the purposes of comparison with aircraft observations and the 

14 use of observed flash rates in the lightning parameterization. Thus, the beginning of the 

15 simulation is assumed to correspond to 1540 UTC. Figure 5 shows a time series of radar 

16 reflectivity plots calculated from the GCE model hydrometeor mixing ratios at 1 km. A 

17 single cell is observed at 30 minutes into the simulation. The cell has completely split in 

18 two at 100 minutes in the model simulation. The southern cell developed a supercell 

19 circulation and has an apparent hook echo at this time. At 150 minutes, a third cell has 

20 developed between the original two cells. The southern cell has begun to decay at 180 

21 minutes, and the northern cell becomes dominant. 

22 Some simulated storm features are similar to observations and some differ. The 

23 model correctly predicted the splitting of the initial cell. However, in the observed 

24 system the northern cell weakened rapidly after the cell splitting event and the southern 
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1 cell became the dominant feature, evolving into a supercell. The simulated southern cell, 

2 while demonstrating supercell characteristics and dominating for a period of time, did not 

3 persist as long as observed. Cloud top heights reached 14 km which compares favorably 

4 with observations [Holler et al., 2000] and the MM5 simulation presented by Fehr et al. 

5 [2004]. Discrepancies between the simulated storm and observations may be because the 

6 nonuniformity of terrain and initial conditions [Stenchikov et al., 2005] were not 

7 accounted for in the GCE simulation. Boundary conditions may also have contributed to 

8 these differences. However, comparison with observations showed that simulated storm 

9 evolution is fairly reasonable for the period of 180 minutes that was chosen for the 

10 chemistry-transport calculations in this study. 

11 At 1657 UTC, just after the cell splitting event, the southern cell was observed by 

12 dual-Doppler radar and the 3-D wind field reconstructed. At this time, a maximum 

13 updraft speed of 24 m s' 1 was observed while the strongest downdraft was 9 m s' 1 (Holler 

14 et al., 2000). At the corresponding time in the simulation (80 minutes), the maximum 

15 updraft velocity was 34 m s' 1 while the maximum downdraft was 7 m s' 1 . Due to the 

16 location of the storm with respect to the radars, dual -Doppler analysis was not possible at 

17 other times. Maximum updraft velocities were approximately 36 m s' 1 between 90 and 

18 130 minutes in the simulation, and then decreased. This is lower than the maximum 

19 updraft of 49 m s' 1 reported in the Fehr et al. [2004] simulation. Downdraft velocities 

20 were also less than those presented in Fehr et al. [2004] in which a maximum of 25 m s' 1 

21 was recorded. Throughout the GCE simulation, downdraft velocities were typically less 

22 than 10 m s' 1 . Low level inflow to the storm occurred between 0.5 and 3 km while 

23 outflow from the anvil was greatest between 9 and 1 1 km. 

24 3.2. Passive Model Results 
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1 The passive version of the CSCTM was used to calculate the transport of CO 2 , 

2 NO x , and O 3 . CO, which in polluted regions has a stronger vertical gradient and would 

3 be preferable to CO 2 as a tracer of upward transport, was not measured on this day. An 

4 initial profile of CO 2 data was constructed using data from the Falcon ascent and a value 

5 of 355 ppbv above the tropopause taken from Strahan et al. [1998]. An initial profile of 

6 O 3 was constructed from data from the Falcon ascent, the DO-228 boundary layer data 

7 for the day, and a climatological average ozone profile for the latitude of Munich above 9 

8 km. The NO x profile was composed of data taken from the Falcon ascent in the free 

9 troposphere and from a profile one standard deviation greater than the average NO x 

10 boundary layer profile during the EULINOX project [Huntrieser et al., 2002]. A profile 

1 1 with values larger than the project mean in the boundary layer was assumed because no 

12 actual measurements were available and high measured boundary layer values of CO 2 

13 and O 3 suggested polluted NO x conditions on this day. Sensitivity calculations were also 

14 performed using the EULINOX boundary layer average and the boundary layer average 

15 plus two standard deviations in order to assess the effects of this uncertainty on the 

16 lightning NO production results. Initial condition profiles for CO 2 , NO x , and O3 are 

17 shown in Figure 6 . IC and CG flash rates for the northern and southern cells (as shown in 

18 Figure 2) were read into the model at 3-minute intervals beginning 21 minutes into the 

19 simulation to correspond to the time when lightning was first observed, and the amount 

20 of lightning NO x produced in each time interval is calculated. 

21 3.2.1. CO 2 and O 3 results 

22 Figure 7 shows a vertical cross section of CO 2 through the southern cell at 150 

23 minutes when the cell was at maximum strength, oriented 65° counterclockwise from due 

24 east. Air containing the maximum CO 2 mixing ratios exceeding 370 ppmv initially in the 
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1 1-2 km region has been transported to over 12 km in the core updraft region, and as high 

2 as 10 km in the anvil, indicating strong upward motion. Both the core and the downwind 

3 anvil regions of the storm are largely comprised of air that resided in the boundary layer 

4 prior to convection, while there is little evidence of entrainment of environmental air with 

5 lower CO 2 mixing ratios. The model also suggests downward transport of smaller mixing 

6 ratios of CO 2 in the 8-1 1 km altitude region behind the storm. 

7 In order to compare simulated tracer transport with mixing ratios observed during 

8 the series of seven anvil penetrations, data collected by the Falcon aircraft were averaged 

9 over approximately 1 1 -second intervals to yield a spatial resolution equivalent to the 

10 model, and then binned into 0.5 km thick layers. Unfortunately, in-cloud observations 

1 1 were only available for three 0.5 km thick layers centered at 8, 8.5 and 9 km AGL. 

12 Therefore, the comparison with model results includes only a 1.5-km thick layer. The 

13 analysis would have benefited from observations at a wide range of altitudes as were 

14 available for the July 12 STERAO-A storm simulated by DeCaria et al. [2005] in which 

15 the aircraft executed a spiral ascent through the storm anvil, measuring NO from 7 to 1 1 

16 km MSL. The area covered by each penetration was calculated from flight position data. 

17 The average distance covered during 6 of the 7 anvil penetrations was determined to be 

18 approximately 24 km in the x-direction and 36 km in the y-direction. A box of this size 

19 was placed around the core of the southern cell where radar observations and flight data 

20 show the Falcon was primarily sampling. The grid cells within this box were sampled at 

21 times in the simulation corresponding to the times of the aircraft sampling at each level, 

22 and cumulative probability distribution functions (pdfs) were calculated for each level. 

23 In addition, the mean, mode and standard deviation of observed and simulated NO x , O3 

24 and CO 2 mixing ratios at 9 km were calculated and are shown in Table 3. The model did 
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1 an excellent job in estimating the mean values at this altitude. However, it appears that 

2 the distribution of observed 03 mixing ratios at this level is substantially broader than 

3 that simulated. 

4 Figure 8 shows the calculated pdfs of observed and simulated CO2 at 8, 8.5 and 9 

5 km AGL. At all three levels, the distribution of simulated CO2 matches the observed 

6 distribution well with a slight overestimation of the maximum values. At 9 km, the 

7 model also underestimates the minimum values, suggesting that the downward transport 

8 at the rear of the storm may not have been as pronounced as seen in the model. Figure 9 

9 shows pdfs of observed and simulated O3 at 8 and 9 km, the only two levels for which a 

10 sufficient number of observations were available to calculate pdfs. At both levels, the 

1 1 simulations underestimate the maximum values and overestimate the minimum values. 

12 The overestimation of minimum values is due in part to the initial condition profiles used. 

13 Though ozone as low as 63 ppbv was observed in the storm, the lowest value in the initial 

14 condition profile was 67 ppbv at 3.5 km because there was no observational evidence 

15 outside the storm to suggest that values lower than this would be appropriate. When 

16 chemical reactions are included in the model (see Section 3.3.) a small loss of ozone 

17 occurs at 8 km, slightly improving the comparison between the simulated and observed 

18 minimum values though overestimation of the minimum values is still noticeable, 

19 particularly at 9 km. The underestimation of the maximum values suggests the model 

20 may be underrepresenting downward transport, although this is not supported by analysis 

21 of the CO2 distributions. This contradiction may be the result of a lack of sufficient 

22 observations of CO2 and O3 near the tropopause to well define the vertical gradients in 

23 this region that are used in the initial condition profiles. Additionally, the observed 

24 discrepancy could be caused by time interpolation of the driving field and approximation 
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1 errors. If downward transport is slightly underrepresented by the model, it is unlikely to 

2 significantly affect the estimates of lightning NO x production because NO x mixing ratios 

3 immediately above the tropopause are similar to the enhanced mixing ratios of NO x 

4 between 7 and 9.5 km in the initial condition profile. 

5 3.2.2. NO x results 

6 To calculate lightning NO x production, Pcg was estimated to be approximately 

7 360 moles of NO per flash based on observed peak current and a relationship between 

8 peak current and energy dissipated from Price et al. [1997]. Several different values of 

9 the Pic/Pcg ratio were simulated and the results compared with observations. The 

10 common assumption that Pic is one tenth Pcg from Price et al. [1997] was simulated and 

1 1 the pdf of observed and simulated NO x at 9 km is shown in Figure 10a. The assumption 

12 that IC flashes are significantly less productive of NO than CG flashes clearly 

13 underestimates NO x at all levels. Fehr et al. [2004] found a Pic/Pcg ratio of 1.4 most 

14 appropriate for a simulation of the same storm. Figure 10b shows the 9 km pdf with NO x 

15 production based on this assumption. At this level, where the majority of NO x 

16 observations were taken, assuming a Pic/Pcg ratio of 1.4 results in an overestimation of 

17 the lightning NO x source. At 8 and 8.5 km, fewer observations are available, but the 

18 maximum observed mixing ratios exceeding 15 and 20 ppbv were reasonably simulated 

19 at these levels using the ratio of 1.4. Therefore the comparison of the Pic/Pcg =1.4 

20 scenario with observations at these two levels is better than at 9 km. 

21 A scenario in which Pic is equal to Pcg was also simulated and the pdfs for 8, 8.5, 

22 and 9 km are shown in Figure 11. At 9 km, the comparison between the observed and 

23 simulated distributions is much better than in the Pic/Pcg of 0.1 and 1.4 scenarios shown 

24 in Figure 10. At 8 and 8.5 km, the model is able to reproduce the distribution below 6 
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1 ppbv fairly well, but fails to produce the large NO x mixing ratios observed. Figure 12a 

2 shows a plot of NO x at 9 km 180 minutes into the Pic=Pcg simulation with the box used 

3 for sampling model output. NO x mixing ratios exceeding 9 ppbv are evident in the core 

4 of the southern cell and mixing ratios over 2 ppbv extend outward in the anvil a distance 

5 of over 70 km. Figure 13 shows the average profiles of NO x mixing ratio within the 

6 averaging box at the beginning of the simulation, after 180 minutes of simulation of NO x 

7 transport without lightning NO production, and at 180 minutes with the inclusion of 

8 lightning NO production. The transport-only simulation results in a mean profile over the 

9 averaging box which maximizes at only 1 ppbv at anvil levels, but this represents an 

10 approximate doubling of the amount of NO x at this level compared with the initial 

1 1 profile. The inclusion of lightning NO x in the model results in an increase in this anvil 

12 maximum to approximately 4 ppbv at 10 km with a second peak of 4 ppbv at 5.5 km. 

13 To determine which production scenario was the most appropriate, the mass of N 

14 in NO x in the column between 7.75 and 9.25 km was also calculated by averaging 

15 observations and model results in each of the three 0.5 km layers. The observations yield 

16 a column mass of 1.13 x 10 3 g N m' 2 . The accuracy of the NO and NO 2 instruments are 

17 5 and 10% respectively [Huntrieser et al., 2002]. Therefore, in terms of measurement 

18 error in the column mass estimate, 10% would be an upper limit. However, there is 

19 additional uncertainty because it is impossible to know how well the aircraft observations 

20 represent a particular area within the storm. Assuming a production scenario in which an 

21 IC flash produces only one tenth as much NO as a CG flash greatly underestimates 

22 column mass, producing 3.43 x 10' 4 g N m’ 2 . The assumption from Fehr et al. [2004] that 

23 an IC flash produces 1.4 times as much NO as a CG flash leads to a column mass of 1.36 

24 x 10" 3 g N m' 2 , an overestimation of approximately 20%. Of the three scenarios 
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1 presented, assuming an IC flash produces as much NO as a CG flash provides the best 

2 comparison with observations with a column mass of 1.05 x 10' 3 g N m' 2 , which 

3 underestimates the column mass calculated from observations by approximately 7%. The 

4 inclusion of chemical reactions in the model tends to decrease NO x , slightly increasing 

5 the error in column mass of the production scenario in which Pic=Pcg=360 moles NO per 

6 flash to 10.5% (see Section 3.3). Based on the comparison of the pdfs and column mass 

7 of the observed and simulated storms, this scenario is selected as the most appropriate of 

8 the three for this storm. An increase or decrease of one standard deviation in the 

9 boundary layer NO x resulted in a change of only 3% in the calculated column mass of N 

10 in NO x . Therefore the assumption of boundary layer NO x mixing ratio has little impact 

11 on our ability to deduce the appropriate lightning NO x production scenario. 

12 A second parameterization scheme was used to estimate NO production per meter 

13 flash channel length. In this parameterization, various values of production per meter are 

14 specified, and production per IC and CG flash is calculated using hourly average 

15 interferometer-observed flash lengths of 21.5, 27.9, and 31.4 km from Thery et al. [2000] 

16 for the hours beginning at 16, 17 and 18 UTC, respectively. By calculating pdfs (Figure 

17 1 1) and column mass, a production of 1.25 x 10" 2 moles NO per meter of flash channel 

18 length was found to yield results comparable to the Pic=Pcg=360 moles NO per flash 

19 scenario, slightly underestimating column mass by approximately the same amount. 

20 Note that the pdfs are nearly identical to those from the first parameterization. 

21 The average hourly flash lengths of Thery et al. [2000] did not differentiate 

22 between IC and CG flashes. Dotzek et al. [2000] attempted to estimate typical flash 

23 lengths for IC and CG flashes separately based on the heights of the main charge layers in 

24 the storm and the diameter of the storm with radar reflectivity greater than 30 dBZ. 
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1 Using this method, they found typical lengths of 43 km for an IC flash, 26.5 km for a 

2 negative CG flash, and 29.5 km for a positive CG flash. It should be noted that because 

3 these estimates are based not on calculated lengths of IC and CG flashes, but on other 

4 parameters, there is a high degree of uncertainty. If these numbers are used in the model 

5 instead of the average hourly flash lengths from Thery et al. [2000], then a production 

6 scenario in which an IC flash produces 8.34 xlO' 3 moles NO per meter of flash channel 

7 length and a CG flash produces 1.35 x 10" 2 moles of NO per meter of flash channel length 

8 is needed to produce a favorable comparison with observations. 

9 3.3. CSCTM (with chemistry) results 

10 To investigate the impact of chemistry on the concentrations of species of interest 

1 1 such as NOx and ozone, a CSCTM run including chemical reactions was performed. 

12 Profiles of CzHe, C2H4, C 3 H6, C 3 H 8 , CH 3 OOH, CO, H 2 0 2 , HCHO, HN0 3 , isoprene, and 

13 PAN were taken from a July mean profile for the appropriate latitude and longitude of the 

14 EULINOX storm computed by the 3-D global University of Maryland Chemical 

15 Transport Model (UMD-CTM) [Park et al., 2004] and are shown in Table 4. The NO to 

16 N0 2 ratio at each CSCTM model level was based on the ratios from the UMD-CTM, but 

17 the initial NO x was equivalent to the values used in the passive version of the model. 

18 Profiles of hydrocarbons were scaled with the aid of airborne hydrocarbon measurements 

19 collected during the 1999 Konvektiver Transport von Spurengasen (KONVEX) campaign 

20 to ensure they represented values typical of the relatively polluted German atmosphere. 

21 Boundary layer concentrations of isoprene were held constant (e.g. at 1 ppbv in the 

22 bottommost layer of the model) during daylight hours to reflect a balance between the 

23 emissions and reactive losses of these compounds. At the conclusion of the 180-minute 

24 simulation, isoprene mixing ratios were approximately 6 pptv in the core updraft region 
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1 of the storm, and were typically less than 1 pptv in the anvil region. A 15-minute “spin- 

2 up” simulation was performed using a column model which included the same chemical 

3 reactions as the full version of the 3-D CSCTM in order to allow the species to come into 

4 equilibrium. The CSCTM was run with the parameterization which uses observed flash 

5 rates as input assuming the lightning NO x production scenario Pic=Pcg=360 moles 

6 NO/flash. Column mass and pdfs were computed from the model output using the same 

7 methods as for the passive version of the model. The use of global model output along 

8 with observed NO x results in a small decrease in NO x in the initial conditions during 

9 spin-up, such that the column mass for the flash rate and flash length scenarios simulated 

10 is 1.01 x 10' 3 g N m' 2 , a difference of 10.5% from the observations. A case in which no 

1 1 lightning NO x is included was also simulated to determine the lightning NO x effects on 

12 in-cloud chemistry. 

13 Figure 12b shows ozone concentrations 180 minutes into the simulation at the 9 

14 km level. Lower ozone air has been transported upwards in the convective cores and is 

15 present in the outflow of the storms. There is significant downward transport of ozone 

16 surrounding the cores of the cells, elevating ozone levels above those of the 

17 environmental air in the anvil region. At the conclusion of the 180 minute simulation, 

18 lightning NO x results in additional O3 production less than 0.1 ppbv in regions outside of 

19 the cloud at 9 km. At 8 and 9 km, the inclusion of chemical reactions in the model 

20 continues to result in the underestimation of maximum ozone mixing ratios (see Figure 

21 9), suggesting that downward transport may be underrepresented in the model. 

22 In order to identify the effect of lightning NO x on O3 production during the storm, 

23 O3 concentrations are averaged at 180 minutes in the simulation within the sampling box 

24 shown in Figure 12a at each model level for both the lightning and no-lightning cases. 
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1 The values from the simulation without lightning NO x are subtracted from the values 

2 from the simulation which included a lightning NO x source. Figure 14 shows the average 

3 change in ozone due to lightning NO x , as well as the maximum and minimum change at 

4 each level. During the lifetime of the storm, the injection of lightning NO x results in a 

5 net loss of ozone averaging less than 4 ppbv at all levels. The maximum net loss during 

6 this 3 hour period exceeds 9 ppbv at 5.5 km. This is due to the large quantities of NO x 

7 (up to 9 ppbv) being introduced into the model. Large NO mixing ratios from lightning 

8 rapidly destroy ozone through the NO + O 3 reaction as described by Wang and Prinn 

9 [2000]. The ozone destruction resulting from including lightning NO x in the model is 

10 likely short-lived. After the cloud dissipates, much of the NO 2 produced by the NO + O 3 

1 1 reaction will be photolyzed to produce NO and 0( 3 P), resulting in O 3 production (see 

12 following section). 

13 3.4. Chemistry-Only Model Results 

14 The chemistry-only version of the CSCTM was used to estimate downstream 

15 ozone production in the 24 hours following the storm. Three-dimensional chemical fields 

16 at 180 minutes in the CSCTM simulation were used to initialize the chemistry-only 

17 version. For these calculations, the storm is assumed to have dissipated and clear-sky 

18 photolysis rates are used. The same 24 x 36 km sampling box shown in Figure 12a was 

19 used to analyze the results at the end of the 24 hour simulations. Table 5 gives the 

20 average mixing ratios within the sampling box at 9 km AGL for a number of species at 

21 the beginning and end of the chemistry-only simulation that included lightning NO x . 

22 Ozone production averaged 1.5 ppbv day 1 in the convective plume at this altitude, while 

23 substantial NO x conversion to HNO 3 took place. Decreases in HCHO and CH 3 OOH 

24 were caused by photolysis. 
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1 The impact of lightning NO x on ozone was examined by averaging the change in 

2 ozone mixing ratios within the box for the lightning and no-lightning cases. Figure 15 

3 shows the change in net ozone production in the 24 hours following the storm due to 

4 lightning NO x , calculated by subtracting the 24-hour change in O 3 from the no-lightning 

5 simulation from the 24-hour change in O 3 from the simulation with lightning NO x . On 

6 average with lightning, there is additional net O 3 production maximizing at 

7 approximately 5 ppbv day ' 1 at 5.5 km. The injection of lightning NO x causes a decrease 

8 in net ozone production averaging less than 2 ppbv day ' 1 between 8 and 10.5 km. There 

9 is a maximum decrease in net O 3 production exceeding 5 ppbv day ' 1 at 9 km, and 

10 maximum additional net production of approximately 9 ppbv day ' 1 at 4.5 km due to 

1 1 lightning . Figure 16 shows a scatter plot of lightning NO x versus the 24-hour change in 

12 net O 3 production resulting from the inclusion of lightning NO in the model for the grid 

13 cells contained in the sampling box at 9 km. The general shape of the plot shows the 

14 change in net ozone production maximizing with lightning NO x mixing ratios less than 1 

15 ppbv, then becoming less positive as lightning NO x increases. After lightning NO x 

16 mixing ratios exceed approximately 2 ppbv, lightning NO x causes decreased net ozone 

17 production in the model. 

18 4. DISCUSSION 

19 The estimated production per CG flash of 360 moles NO from the CSCTM for 

20 this storm is close to the value of 330 moles NO per flash from Fehr et al. [2004], with 

21 the difference due to the criteria for designating which BLIDS flashes were considered 

22 actual CG flashes in the two models. However, the assumed Pic/Pcg ratio of 1.4 yields a 

23 production per IC flash of 462 moles NO which overestimates NO x when compared to 

24 the observed pdfs and column mass of observations. The production per CG flash of 360 
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1 moles NO is less than the production estimated for the July 12 STERAO storm by 

2 DeCaria et al. [2005] of 460 moles NO per CG flash. The IC production rate of 360 

3 moles NO per flash is in the range of 345 to 460 moles NO per IC flash estimated by 

4 DeCaria et al. [2005] on the basis of calculated column mass and the average vertical 

5 NO x profile. Assuming a global average flash rate of 44 flashes s' 1 [Christian et al., 

6 2003], the lightning NO x production scenario from the EULINOX storm yields an annual 

7 global lightning NO x source of 7 Tg N yr' 1 . Huntrieser et al. [2002] also estimated the 

8 global lightning NO x production rate in two ways. The first method involved calculating 

9 the average anvil NO x mixing ratios during EULINOX storms and multiplying by the 

10 average air flux out of the anvil and an estimate of the average number of thunderstorms 

1 1 occurring globally. The second method involved estimating NO production per meter 

12 flash channel length from NO spikes observed in the July 21 EULINOX storm and 

13 assuming a mean flash length and global flash rate. These methods yielded estimated 

14 global production rates of 3 and 4 Tg N yr' 1 , respectively both of which are less than the 

15 estimate of 7 Tg N yr' 1 presented above. These differences could be because other 

16 EULINOX storms may have been less productive of lightning NO than the July 21 storm, 

17 or because the NO spikes analyzed from the July 21 storm may not have been 

18 representative of the lightning NO production in the storm. Our estimate of 1.25 x 10' 2 

19 moles NO per meter flash channel length also compares favorably with other estimates of 

20 NO production per meter available in the literature. This estimate is within the range of 

21 3.3 x 1 O' 4 and 1.7 x 10' 2 moles NO per meter flash channel length from Stith et al. [1999] 

22 based on the July 10 STERAO storm. It is significantly larger than the value of 1.7 x 10' 3 

23 moles NO per meter flash channel length estimated by Skamarock et al. [2003] which 

24 included a number of short duration interferometer flashes not included in other studies 
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1 of STERAO cases (e.g. DeCaria et al, 2005). It is also larger than the estimate of 4.5 x 

2 10' 3 moles NO per meter flash length obtained by Huntrieser et al. [2002]. 

3 DeCaria et al. [2005] also examined the change in the net production of ozone 

4 resulting from lightning NO x in the convective plume following convection and found an 

5 average anvil net ozone production increase maximizing at 10 ppbv day' 1 at 9 km AGL. 

6 Maximum increased ozone production at this level exceeded 12 ppbv day' 1 . The 

7 simulations presented of the July 21 EULINOX storm show on average smaller changes 

8 in net ozone production, maximizing at 5 ppbv day' 1 at 5.5 km AGL. The maximum 

9 increase in the net production of ozone was greater than 9 ppbv day' 1 at 4.5 km. The 

10 EULINOX storm contained much larger NO x mixing ratios than the STERAO storm due 

11 to convective transport of boundary layer pollution as well as more lightning. Assuming 

12 a Pic/Pcg ratio of 1 for the STERAO storm yielded an average in-cloud NO x mixing ratio 

13 of only 0.9 ppbv at 9 km AGL (DeCaria et al., 2005), while the average at 9 km AGL 180 

14 minutes into the simulation is 3.3 ppbv in the EULINOX storm. As a result, NO x was 

15 less efficient at O3 production in the EULINOX storm. The decrease in ozone production 

16 seen in the simulation of the EULINOX storm itself may be accentuated due to the 

17 changes to the lightning NO x parameterization which attempt to simulate individual 

18 lightning flashes and the spiky nature of NO x observations by placing NO from a flash 

19 into a small subset of randomly selected grid cells rather than uniformly distributing it 

20 over a larger area of the simulated cloud (as implemented by DeCaria et al., 2005). The 

21 higher NO x mixing ratios cause a loss of ozone at some levels and only minimal amounts 

22 of production of ozone at others within the storm. 

23 Ozone production following convection has also been studied in global chemical 

24 transport models. Park et al. [2004] used a global stretched-grid chemical transport 
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1 model (UMD-CTM) which included parameterized convection and lightning NO x 

2 production to examine the change in ozone following a convective event observed over 

3 Kansas and Oklahoma in 1985. When comparing ozone mixing ratios prior to convection 

4 with ozone mixing ratios in the day following convection, they found an enhancement in 

5 the downstream region at 9 km averaging 7 ppbv and ranging from zero to 25 ppbv. 

6 Prior to convection, O 3 at 9 km in the EULINOX simulation was 80.1 ppbv and after the 

7 conclusion of the 24-hour chemistry-only simulation, O 3 mixing ratios within the 

8 sampling box ranged from 80.9 to 97.2 ppbv (an increase of ~1 to 17 ppbv) with an 

9 average of 89.1 ppbv, an increase of 9 ppbv. The range of the increase in O 3 following 

10 convection calculated in the EULINOX simulation is less than the range of enhancement 

1 1 from Park et al. [2004]. The average increase in O 3 of 9 ppbv is slightly greater than the 

12 7 ppbv average from the Park et al. [2004] study. 

13 5. SUMMARY 

14 The 3-D GCE cloud resolving model initialized with a single sounding was used 

15 to simulate a powerful thunderstorm observed on July 21, 1998 as part of the EULINOX 

16 field campaign. The model successfully reproduced a number of observed features of the 

17 storm including the splitting of the storm into two distinct cells, one of which developed 

18 supercell characteristics, and the other of which became multicellular. Time series of IC 

19 and CG flash rates were generated for both cells using BLIDS and interferometer data, 

20 and used in a parameterization of lightning NO production included in an offline 3-D 

21 CSCTM simulation driven by fields from the GCE simulation. This newly developed 

22 parameterization provides a means to better simulate the NO x resulting from lightning 

23 flashes in the highly electrified portion of a storm than is possible with the DeCaria et al. 

24 [2005] scheme. 
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1 Probability distribution functions were calculated for CO 2 , O3, and NO x from 

2 aircraft observations taken during 7 penetrations of the storm and model output. A 

3 comparison of CO 2 and O 3 pdfs showed the model adequately represented the transport 

4 of tracer species. However, the lack of sufficient observations in the upper troposphere 

5 and lower stratosphere for use as initial conditions hampered a full evaluation of the 

6 model transport. Several lightning NO x production scenarios were simulated and the 

7 results were compared with observations to determine which was most appropriate for the 

8 storm. The scenario in which IC and CG flashes both produce 360 moles of NO per flash 

9 compared most favorably with observations, both in terms of pdfs and column mass. 

10 This production scenario for a strong supercell over Germany is equivalent to an annual 

1 1 global source of 7 Tg N yr' 1 . A second lightning NO x parameterization using average 

12 hourly flash lengths was employed to estimate a production per meter flash channel 

13 length of 1.25 x 10 2 moles NO. 

14 Many modeling studies have assumed IC flashes are significantly less productive 

15 of NO than CG flashes because they are less energetic and pass through lower pressure 

16 layers. Price et al. [1997] proposed that an IC flash typically produces 10% as much NO 

17 as a CG flash based on estimates of energy dissipation. This has been disputed by Zhang 

18 et al. [2003] who estimate that IC flashes dissipate between 50 and 100% of the energy of 

19 CG flashes using recent observational evidence and a reanalysis of the assumptions made 

20 in the Price et al. [1997] study. Though a typical IC flash might be expected to produce 

21 somewhat less NO than a typical CG flash, the differences in pressure where the two 

22 types of flashes typically occur and in energy dissipation may be effectively countered by 

23 the fact that IC flashes can have greater lengths. Dotzek et al. [2000] estimated that the 

24 average length of an IC flash was greater than the average length of a CG flash in the July 
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1 21 EULINOX supercell based on lightning and radar observations. Though IC and CG 

2 flashes on average may produce equivalent quantities of NO, our results suggest that an 

3 IC flash may be less productive of NO per meter flash channel length than a CG flash. 

4 Assuming either the production scenario per flash or per meter flash channel 

5 length leads to a loss of ozone at all levels during the lifetime of the storm due to the 

6 resulting large NO mixing ratios. Chemical fields from 180 minutes in the CSCTM 

7 simulation were used to initialize a chemistry-only version of the model which includes 

8 diffusion but no further cloud-scale advection. The entire model domain is assumed to be 

9 translated downstream by the wind. This model was integrated forward 24 hours to 

10 estimate the effect of lightning NO x production on ozone in the day following the storm. 

1 1 Layer average increases in net production of ozone were found to maximize at 5 ppbv 

12 day' 1 at 5.5 km, though grid cells (particularly at 8 to 10.5 km) with the highest 

13 concentration of lightning NO resulted in decreased net ozone production. We speculate 

14 that outflow from highly electrified storms or storms occurring over heavily polluted 

15 regions in other parts of the world (e.g. eastern U.S., east Asia, etc.) may also contain air 

16 parcels with reduced O3 production. Dilution of these convective plumes as they are 

17 advected farther downwind will result in a transition to increased O3 production. It 

18 appears that net O3 production over the first 24 hours in outflow from storms with large 

19 flash rates over polluted regions may be less than in the outflow from moderate flash rate 

20 storms over cleaner regions. 
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20 Figure Captions 

21 Figure 1. 1° elevation scans from the POLDIRAD radar at 1640, 1652, 1733, and 1802 

22 UTC with BLIDS CG flashes recorded during the 2 minutes before and after each radar 

23 scan overlaid. Horizontal arrows denote positive flashes and jagged arrows indicate 

24 negative flashes [Holler et al., 2000]. 
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1 


2 Figure 2. Time series of IC and CG flash rates for the (a) northern and (b) southern cells. 

3 

4 Figure 3. Locations of VHF sources recorded by interferometer at stations ST1 and ST2 

5 between 1740 and 1810 UTC overlaid with the Falcon’s trajectory during the same time 

6 period. Red letters indicate the location of electrified cells [Huntrieser et al., 2002]. 

7 

8 Figure 4. (a) Vertical distribution of the number of grid cells in the horizontal included at 

9 each model level in a CG flash (solid line) and in an IC flash (dotted line), (b) Schematic 

10 diagram of the lightning parameterization. Dashed line represents computed 20 dBZ 

1 1 radar reflectivity contour at 9 km 150 minutes into the simulation. Triangle shows the 

12 location of the maximum updraft and the larger box is the area from which an initiation 

13 point for the lightning is selected. The circle marks the initiation point, and the smaller 

14 box shows the area in which the flash is constrained. 

15 

16 Figure 5. Radar reflectivity at 1 km elevation computed from GCE hydrometeor fields at 

17 (a) 30, (b) 100, (c) 150, and (d) 180 minutes in the simulation. 

18 

19 Figure 6. Initial condition profiles of (a) CO2, (b) NO x , and (c) O3. 

20 

21 Figure 7. Cross section of CO2 mixing ratios in the southern cell from the CSCTM at 150 

22 minutes at an angle of 65° counterclockwise from east. 

23 
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1 Figure 8 . Pdfs of observed (solid) and simulated (dashed) C0 2 mixing ratios at (a) 8 , (b) 

2 8.5, and (c) 9 km. 

3 

4 Figure 9. Pdfs of observed (solid) and simulated (without chemical reactions - dashed, 

5 with chemical reactions - dotted) O 3 mixing ratios at (a) 8 and (b) 9 km. 

6 

7 Figure 10. Pdfs of observed (solid) and simulated (dashed) NO x mixing ratios at 9 km 

8 assuming (a) Pic/Pcg = 0 . 1 , and (b) Pic/Pcg = 1.4. 

9 

10 Figure 1 1. Pdfs of observed (solid) and simulated (assuming Pic=Pcg=360 moles NO - 

1 1 dashed, assuming P= 1.25x1 0‘ 2 moles NO per meter flash channel length - dotted ) NO x 

12 mixing ratios at (a) 8, (b) 8.5, and (c) 9 km assuming Pic/Pcg = 1.0. 

13 

14 Figure 12. (a) NO x mixing ratios at 9 km elevation assuming Pic=Pcg=360 moles 

15 NO/flash at 180 minutes in the passive CSCTM simulation. The box indicates the grid 

16 cells sampled for calculation of column mass and pdfs and the black line indicates the 20 

17 dBZ contour of computed radar reflectivity, (b) O 3 mixing ratios at 9 km at 180 minutes 

18 in the CSCTM simulation including chemical reactions. The black line indicates the 20 

19 dBZ contour of computed radar reflectivity. 

20 

21 Figure 13. Mean NO x mixing ratios within the averaging box at the beginning of the 

22 passive model simulation (solid), and after 180 minutes with only transport simulated 

23 (dashed), and with lightning and transport simulated (dotted). 

24 
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1 Figure 14. Change in O 3 mixing ratios due to lightning N0 X during the lifetime of the 

2 storm. Solid line is the average (over the sampling box) and brackets indicate minimum 

3 and maximum change. 

4 

5 Figure 15. Change in net O3 production due to lightning NO x in the 24 hours following 

6 the storm. Solid line is the average (over the sampling box) and brackets indicate 

7 minimum and maximum change. 

8 

9 Figure 16. Change in net O3 production due to lightning NO x in the 24 hours following 
10 the storm versus lightning NO x at 9 km for grid cells contained in the sampling box. 
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1 Table 1. Estimates of NO Production per Flash 


Author(s) 

Field 

Project 

Method 

NO production 
(moles NO per flash) 

DeCaria et al. [2000] 

STERAO 

2-D cloud model, 

aircraft 

observations 

230 - 460 

DeCaria et al. [2005] 

STERAO 

3-D cloud model, 

aircraft 

observations 

345 - 460 

Fehretal. [2004] 

EULINOX 

3-D cloud model, 

aircraft 

observations 

330-462 

Huntrieser et al. [1998] 

LINOX 

Aircraft 

observations 

66 - 498 

Skamarock et al. [2003] 

STERAO 

3-D cloud model, 

aircraft 

observations, 

lightning 

observations 

43 

Thery et al. [2000] 

EULINOX 

Lightning 

observations 

500 
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1 Table 2. Estimates of NO Production per meter Flash Channel Length 


Author(s) 

Field 

Project 

Method 

NO production 
(moles NO per m flash 
length) 

Holler etal. [1999] 

LINOX 

Lightning 

observations 

1.7 - 6.6 x 10' z 

Huntrieser et al. [2002] 

EULINOX 

Aircraft 

observations 

4.5 x 10' 3 

Skamarock et al [2003] 

STERAO 

3-D cloud model, 

aircraft 

observations, 

lightning 

observations 

1.7 x 10‘ 3 

Stith et al. [1999] 

STERAO 

Aicraft 

observations 

3.3 x 10' 4 - 1.7 x 10' 2 
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1 Table 3. Statistics of Observed and Simulated Tracer Mixing Ratios at 9 km 



Mean 

Mode 

Standard Deviation 

Observed CO 2 (ppmv) 

367.0 

366.4 

1.5 

Simulated CO 2 (ppmv) 

367.8 

371.8 

3.4 

Observed O 3 (ppbv) 

90.0 

84.8 

15.8 

Simulated O 3 (ppbv) 

89.1 

84.1 

5.6 

Observed NO x (ppbv) 

2.4 

1.3 

1.7 

Simulated NO x (ppbv) 

2.62 

0.3 

2.1 


2 

3 * Statistics are computed over the sampling box shown in Figure 12a. 


43 




1 Table 4. Initial Condition Profiles of Species Taken from the UMD-CTM and Included 

2 in the CSCTM 


Altitude 

(km) 

C 2 H 6 

(pptv) 

C 2 H 4 

(pptv) 

C 3 H 6 

(pptv) 

c 3 h 8 

(PPtv) 

HCHO 

(pptv) 

hno 3 

(PPtv) 

H 2 0 2 

(pptv) 

PAN 

(pptv) 

1 

1260 

140.0 

200.0 

298.0 

1030.0 

1800 

1280 

776 

2 

820 

32.9 

60.7 

114.0 

452.0 

593 

1050 

352 

3 

718 

15.9 

38.9 

71.4 

308.0 

323 

845 

366 

4 

706 

12.2 

27.1 

63.0 

223.0 

359 

685 

488 

5 

716 

11.3 

25.5 

65.7 

175.0 

314 

612 

587 

6 

736 

13.1 

35.2 

73.5 

151.0 

433 

545 

648 

7 

756 

15.9 

50.2 

81.9 

136.0 

530 

469 

675 

8 

778 

19.0 

65.5 

88.8 

129.0 

593 

409 

680 

9 

804 

22.2 

75.8 

95.4 

132.0 

641 

403 

670 

10 

820 

21.7 

65.4 

98.7 

115.0 

654 

323 

641 

11 

830 

17.8 

45.6 

98.4 

88.1 

649 

229 

606 

12 

826 

7.1 

14.6 

91.2 

49.4 

615 

146 

562 

13 

856 

3.8 

4.3 

96.0 

36.5 

647 

119 

531 

14 

902 

3.6 

2.7 

106.0 

34.4 

713 

118 

508 

15 

958 

4.6 

4.1 

119.0 

36.8 

838 

134 

501 
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1 


Table 5. Average Species Concentrations at 9 km at the Beginning and End of the 
Chemistry-Only Simulation 


Species 

Average Mixing Ratio Immediately 
Following Convection (ppbv) 

Average Mixing Ratio 24 hours 
after Convection (ppbv) 

NO x 

3.28 

2.11 

o 3 

87.6 

89.1 

hno 3 

0.334 

1.25 

HCHO 

0.386 

6.48 x 10 2 

h 2 o 2 

0.524 

0.339 

ch 3 ooh 

0.160 

8.21 x 10' 2 

ch 3 co 3 no 2 

0.650 

0.671 

OH 

1.57 x 10 3 

1.07 x 10' 5 

ho 2 

1.53 x 10' 4 

2.80 x 10‘ 4 

ro 2 

4.57 x 10' 3 

8.77 x 10 5 1 


4 * Averages are computed over the sampling box shown in Figure 12a. 
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1 Figure 2. 

2 
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a) ^ July 21 Flashrates - Northern Cell 
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1 Figure 5. 
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a) Radar Reflectivity(dBZ) Z = 1km 30mln b) Radar Reflectivity(dBZ) Z = 1km lOOmln 
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C0 2 (ppmv) angle = 65 deg. 150 min. 
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Figure 12. 


a) N0 X (ppbv) Z = 9 km 180 minutes 
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